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ABSTRACT 

A five square arcminute region around the luminous radio-loud quasar SDSS 
J0836-I-0054 {z — 5.8) hosts a wealth of associated galaxies, characterized by 
very red (1.3 < ■1775 — ^sso < 2.0) color. The surface density of these 2; ~ 5.8 
candidates is approximately six times higher than the number expected from 
deep ACS fields. This is one of the highest galaxy overdensities at high redshifts, 
which may develop into a group or cluster. We also find evidence for a sub- 
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structure associated with one of the candidates. It has two very faint companion 
objects within two arcseconds, which are hkely to merge. The finding supports 
the resuhs of a recent simulation that luminous quasars at high redshifts lie on 
the most prominent dark-matter filaments and are surrounded by many fainter 
galaxies. The quasar activity from these regions may signal the buildup of a 
massive system. 

Subject headings: quasars: individual (SDSS J0836-I-0054), large-scale structure 
of Universe, galaxies: high-redshift 

1. INTRODUCTION 

The fluctuations in the cosmic microwave background temperature observed by WMAP 
(Bennett et al. 2003) are beheved to be the seeds for the first generation of baryonic objects, 
which eventually evolved into the stars, galaxies and clusters seen in the current epoch. 
Several hundreds galaxy candidates at ^ ^ 6 have been found as i-dropout objects (Stanway 
et al. 2003; Yan & Windhorst 2004; Bouwcns ct al. 2003, 2005), characterized by a large color 
difference between the 2775- and 2;85o-band and faint magnitudes of approximately 26. The 
current best estimate of the average surface density of these galaxies is from the GOODS 
survey (Giavalisco et al. 2004); approximately 0.25 galaxy per square arcminute to a -Zgso- 
band magnitude of 26.5 (Dickinson et al. 2004; Bouwens et al. 2005). In addition, at least 
a dozen of quasars have also been discovered at redshift z ^ 5.7 (Fan et al. 2004), and their 
high luminosities suggest the presence of massive black holes of > 10^ Mq. The formation 
of such black holes at epochs less than one biUion years after the Big Bang requires an 
extremely high accretion rate (Haiman & Loeb 2001). 

According to the hierarchical model of galaxy formation (Press & Schechter 1974), 
massive baryonic objects are formed in the densest halos. The latest simulations of CDM 
growth (Springel et al. 2005) predict that 2; ~ 6 quasars may lie in the center of very massive 
dark matter halos of ~ 4 x 10^^ Mq, which will grow into the most massive clusters of 
~ 4 X 10^^ Mq at z = 0. Studying the clustering properties of young galaxies and quasars 
may provide insights into their formation history. Reports of clustering of galaxies at redshift 
z ~ 6 (Ouchi et al. 2005; Wang et al. 2005; Malhotra et al. 2005) suggest that the formation 
of large structures in the early Universe may be more significant than the cold dark matter 
(CDM) models have predicted. Recently Stiavelli et al. (2005) reported that one of the most 
distant quasars, SDSS J1030-I-0524 at ^ = 6.28, exhibits an excess of associated galaxies. 
It has also been known that some radio galaxies harbor an enhanced number of associated 
galaxies. Venemans et al. (2002), and Miley et al. (2004) found a potential protocluster 
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around TN1338-1942, a radio galaxy ect z — 4.1. Venemans et al. (2004) and Overzier et 

al. (2005) observed an excess of Lya emitters (LAEs) and Lyman break galaxies (LBGs) in 
the field of TN0924-2201, the most distant radio galaxy at z = 5.2. These reports confirm 
an empirical result that powerful, high-redshift radio galaxies are associated with massive 
forming galaxies (van Breugel et al. 1999). To explore the possibility that radio-loud quasars 
may also be a signpost of galaxy clustering at high redshifts, we initiated a project with 
HST/ACS to image the fields around some of the most distant radio-loud quasars. 

The quasar SDSS J0836+0054 {z = 5.82; Fan ct al. 2001) is the most distant radio-loud 
quasar known to-datc, and is one of the most luminous. The object is detected in the FIRST 
radio survey, with a total flux density at 1.4 GHz of 1.11 niJy. VIA and MPMBA-IRAM 
observations (Petric et al. 2003) yield flux densities of 1.75 mJy at 1.4 GHz, 0.58 mJy at 5 
GHz, and a non-detection at 250 GHz, with a 3a upper limit of 2.9 mJy. It is a compact 
steep-spectrum radio source with a radio spectral index of —0.8. VLBI observations at ~ 10 
milliarcsecond angular resolution (Prey et al. 2003) show an apparent core-jet morphology, 
with no indication of multiple images produced by gravitational lensing. The quasar's enor- 
mous power at M ~ —27.8 in the rest-frame UV band implies a mass of the central black 
hole of 5 X 10^ Mq, posing a challenge to the theoretical models on how such massive objects 
are formed in the very early Universe. 

Throughout the paper, we use AB magnitudes and assume the common values of cos- 
mological parameters: (0^,0^) = (0.3,0.7) and Hq = 70 kms~^Mpc~^. 

2. DATA 

The HST/ACS observations of SDSS J0836+0054 were carried out on 2004 October 
8 and November 17, with a total exposure of 10778 sec in the ^850-band and 4676 sec in 
the Z775-band, respectively, at two different position angles. They were processed with the 
standard pipchne CALACS (Pavlovsky ct al. 2005), then with APSIS (Blakeslcc ct al. 2003). 
These procedures carried out flat-fielding, removing bias, dark-current and cosmic-ray events, 
corrections for geometrical distortion of the detectors, and drizzling the dithered images. The 
final images cover approximately 11.4 square arcminutes. The APSIS tasks determine the 
Galactic extinction of E(B-V)=0.05 from the dust maps of Schlegel et al. (1998) and apphed 
0.1 and 0.07 to 2775- and ^gso-band magnitudes, respectively. 

We used SExtractor (Bertin & Arnouts 1996) to find and parameterize sources from the 
science images and their RMS counterparts. We first used the 2;850-band as the detection 
image and then reran the task in a dual mode, namely to use the profile information in 
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the 2;85o-band to link and constrain the parameter counterparts in the i775-band image. The 
limiting magnitudes are similar to those of the GOODS fields: 26.5 in the Zgso-band, for a 
10a detection of a source of &!2 in diameter. As shown in Table 1, we selected objects with 
a color 2775 ~ -2850 > 1-3 (MAGJSO). Only sources with a star-galaxy index of < 0.8 were 
considered (0 for galaxy and 1 for star). To further avoid contaminations from cosmic ray 
events, we only considered sources in the sky region that is covered by all six exposures in 
each band, which is approximately 10 square arcminutes. The limiting zgso magnitudes are 
MAG_AUTO as they allow us to collect most of the source flux, but the colors are determined 
with MAG_ISO, in order to maximize the signal-to- noise (S/N) ratio. Extensive tests by 
Benitez et al. (2004) suggest that the colors of faint galaxies are more robust with MAG_ISO. 
The color selection threshold of 1.3 is chosen to avoid the contamination from low-rcdshift 
interlopers. Seven sources are listed, and five of them are 2775-faint (i-faint) objects, enabling 
us to separate them from objects at z > 6. Lya emission at 6 < 2; < 7 is at wavelengths 
redward of the i775-band, and the 2775 — ^§50 color is at least 2.5 magnitudes. No detection 
in the i775-band is therefore anticipated. There are nine objects with 1.0 < ^775 — 2:850 < 1-3, 
which may be LAEs at 2; ~ 5.8, or galaxies at 2; ~ 1. 

At a redshift of 5.8, one arcsecond corresponds to approximately 5.7 kpc, and the 
transverse dimension of the ACS field is approximately 1.1 Mpc. Along the line of sight, the 
proper distance between z — 5.7 and 5.9 is approximately 13 Mpc. 

3. RESULTS 

We identified seven galaxy candidates with z > 5.5 in the field of this quasar. Table 
1 lists the properties of the candidates, along with that of the quasar, and Fig. 1 displays 
their positions. Object D is near a foreground cluster approximately one arcminute south 
of the quasar. The cutout images of each candidate are shown in Fig. 2. We carried out a 
test on negative images (Dickinson et al. 2004) as spurious sources may be present even at 
S/N > 5. The 2775- and 2;85o-band images were multiplied by —1, and SExtractor tasks were 
carried out using the same detection parameters. No candidate was detected with the same 
selection criteria. 

The Bayesian photometric redshifts (BPZ; Benitez 2000) were calculated with a cali- 
brated template set of Benitez et al. (2004) supplemented by a very blue starburst template. 
These templates significantly improve the photometric redshift estimation for faint, high- 2; 
galaxies. The new features in BPZ identify multiple peaks of redshifts and assigns a proba- 
bility to each of them. The values of the first-peak redshifts are listed in column 8 of Table 
1. 
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While it is common to refer to objects with large ^775 — 2:850 color as "i-dropouts" , we can 
obtain a better redshift discrimination by selecting objects with strong 2775 — zg^o breaks, but 
which are still detected in the i775-band. At 2; ~ 5.8, the redshifted Lya feature straddles 
both the 1775- and ZgsQ-band. Because of the high throughput of ACS, most candidates at 
this redshift are expected to be detected in the i775-band. As the redshift increases from 5.8 
to 6.0, Lya emission rapidly moves out of the i775-band, and therefore the ^775 — zs5o color 
increases rapidly from magnitude ~ 1.5 to ^2.0. Having set a detection limit of Zs5o < 26.5, 
we are able to secure a 2a detection of sources with ^775 — zgso < 2.0. We tested the reality 
of i775-band detections by deliberately shifting the i775-band image by ±20 pixels along its 
row and column and then reran SExtractor. The source apertures defined by the ^gso-band 
images would point to a nearby sky field, and in all the cases, the aperture source fluxes in 
the i775-band were below S/N=1.4. We therefore conclude that an i775-band detection at S/N 
> 2 is unlikely the result of random fluctuations. The detection in both bands significantly 
enhances the confidence level, and also enables us to exclude background galaxies at z > 6 
(4775- or i-dropouts). Objects B and G in Table 1 are not detected in the ^775 band. They 
are likely background galaxies at z > 6 and not physically associated with the quasar. 



4. DISCUSSION 

The GOODS results (Dickinson et al. 2004; Bouwens et al. 2005) suggest that the surface 
density of i-dropouts and i-faint objects is approximately 0.25 per square arcminute, to a 
limiting magnitude of ^850 ~ 26.5. These objects are believed to be at redshifts 5.5 < z < 6.5. 
Since we only select i-faint objects, the photometric redshifts of these objects only cover the 
range of 5.5 < z < 6. Prom previous studies (Dickinson et al. 2004; StiaveUi et al. 2005; 
Bouwens et al. 2005) we estimate that approximately 60% of these objects are i-faint and the 
rest i-dropouts. The surface density of i-faint objects is only ~ 0.15 per square arcminute, or 
1-2 in the ACS field of view. Since the chance of finding seven i-dropouts (including i-faint 
objects) in GOODS in a random ACS/WFC-size cell (~ 11 square arcminutes) amounts to 
a few percent (Stiavclli ct al. 2005), our finding of five i-faint objects (and two i-dropouts) 
in one ACS field suggests a significant source overdensity. Since all our candidates lie in a 
region of five square arcminutes the actual factor of overdensity is approximately six with 
respect to GOODS, with no random cells drawn from GOODS containing seven objects (four 
being the highest using the sample of Bouwens et al. 2005). Although cosmic variance is 
expected to be non-neghgible even on scales as large as probed by GOODS (Somerville et al. 
2004), field-to- field variations can not be determined empirically until larger surveys become 
available. 
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Our results provide new evidence for an excess of galaxies associated with quasars at 
z ^ 5.8. Stiavelli et al. (2005) find seven candidates at ^ ^ 6 in the field of a radio-quiet 
quasar at z = 6.28. Four out of these seven candidates are detected in the i775-band, and are 
likely to be foreground galaxies at redshift 5.5 < z < 6. The number of galaxies associated 
with the z = 6.28 quasar is at most three. We therefore believe that the five z ~ 5.8 
candidates in the vicinity of quasar SDSS J0836+0054 represents a significant overdensity. 
Such an overdensity is consistent with the prediction of the Millennium Simulation (Springel 
et al. 2005) that a 'first quasar' candidate at z — 6.2 lies on one of the most prominent 
dark-matter filaments and is surrounded by a large number of other, much fainter galaxies. 
The quasar itself exhibits an ^775 — 2:850 color of 1.2, which is considerably smaller than 
that in the SDSS (2.2). This is because the i775-band sensitivity of ACS detectors extends 
further to the red than that of the SDSS. The bulk of rcdshiftcd Lya emission falls into that 
extended wavelength region beyond 8000 A and boosts the i775-band flux. We carried out 
color simulations of galaxies at 2; ~ 5.8, using the spectral evolution library of Bruzual & 
Chariot (2003). Among the 39 representative model templates, we choose 21 at age of 1.4 
Gyr and younger, which include star-burst models. We also added Lya emission lines to the 
models with ages of 25 Myr. The Lya emission is assumed to have a hue width FWHM=1000 
kms"^, and varying equivalent widths 100 or 200 A in the rest frame. All the model spectra 
are redshifted, and corrected for intergalactic absorption. As shown in Fig. 3, the ^775 — zgso 
color peaks around 1.7 for LBGs, and that of LAEs extends at 1.0 < ^775 — 2850 < 1-6. 
Approximately 90% of the LBGs are selected by our criteria, but only a fraction for LAEs. 

The field around source C is interesting. The source seems to be a multiple system, as 
shown in Fig. 4. In its vicinity there are two red objects, which are slightly fainter and are 

marked as C2 and C3. The two companion objects of source C are not among our initial 
sample, as they do not meet our selection criteria. However, they are in the vicinity of 
source C and share similarly red colors suggesting that they may be physically associated. 
The angular separation between these three sources is only ~ 1.7 arcsecond or ~ 10 kpc, 
suggestive of merging. 

There is a slight chance that some of the candidates are reddened objects at lower 
redshifts. Red galaxies at 2 ~ 1 display a color of 0.9 ^ ^775 — 2850 ^ 1-1 (Mei et al. 2005). 
Uncertainties in the ^775 — 2850 color in our sample may introduce contaminations from red 
galaxies at 2; ~ 1. While it is possible that large uncertainties in magnitudes may lead to 
contamination from these galaxies, such an effect should also apply to the ACS images of 
the GOODS fields because of the same instrument and comparable exposure times. As the 
GOODS fields exhibit a considerably lower surface density of i-faint objects than our field, 
the contamination level is low at ~ 8% (Bouwens et al. 2003). The so-called Extremely 
Red Objects (ERO) can display sharply rising fluxes toward longer wavelengths in several 
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adjacent bands. The lack of deep infrared images does not allow us to exclude such objects, 
but since the sizes and morphologies of our candidates do not fit either type, having neither 
pointlike or rather diffuse profiles, we conclude that they are unlikely EROs. Malhotra et al. 
(2005) suggest that ~ 90% of sources with ^775 — -^850 > 1-3 and 2:859 < 27 are spectroscopically 
confirmed at redshift z > 5.5. The actual contamination rate is even lower when star-like 
objects are excluded. 

Our finding complements a number of recent reports of overdensities of associated galax- 
ies in the vicinities of radio galaxies at 2 > 4 (although not exclusively) and suggest a possi- 
bility of enhanced activities of clustering in the fields of radio sources at high rcdshifts. The 
galaxy candidates are more than a hundred times fainter than the quasar itself, and they 
form a good example of a hierarchical evolution. The detection of radio emission signals 
an environment that is rich in dark matter, thus harboring the formation of massive black 
holes. 



5. SUMMARY 

The region around the radio-loud quasar SDSS J0836-I-0054 is rich in i-faint objects, 
which are candidates for galaxies at z ~ 5.8. The surface density of these objects is ap- 
proximately 4-6 times higher than that of the GOODS fields, yielding one of the highest 
overdensities at 2; ~ 6 known to-date. Our finding supports a hierarchical structure pre- 
dicted by a recent simulation that luminous quasars arc surrounded by fainter galaxies. 
Spectroscopic observations are needed to further confirm the association and enable us to 
estimate the volume density occupied by these sources. Future observations of this field in 
infrared bands will also provide information about the spectral energy distribution of these 
distant sources. 

The observations add fresh evidence that quasars and radio galaxies are good beacons for 
finding protoclusters of young galaxies at high rcdshifts. Our measurements could provide 
the first constraint on the halo occupation number for LBGs in one of the most massive 
haloes at high redshift, which will provide interesting comparison to numerical simulations. 

ACS was developed under NASA contract NAS 5-32865, and this research has been 
supported by NASA grant NAG5-7697 and by an equipment grant from Sun Microsystems, 
Inc. The Space Telescope Science Institute is operated by AURA Inc., under NASA contract 
NAS5-26555. We are grateful to K. Anderson, J. McCann, S. Busching, A. Pramarini, 
S. Barkhouser, and T. Allen for their invaluable contributions to the ACS project at JHU. 
We thank the anonymous referee for constructive comments. 
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Table 1. Objects with Large ^775 — ^850 Color 
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^Calculated using FLUXJSO and MAGJSO. 

''Bayesian photometric redshift, at a 67% confidence level. First-peak redshifts are estimated with a 
preset upper limit oi z — 7.0. 

'^i-dropout. Not considered as being associated with the quasar. 

"^Calculated using MAG_AUTO, as MAG_ISO yields no detection. 
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Figure Captions 

Fig. 1.— Composite HST/ACS image of the field of the quasar SDSS J0836+0054. The field 
size is approximately 3 arcminutes (~ 1.1 Mpc of proper distance). The quasar (Q) and the 
candidates are marked with open circles. 

Fig. 2. — Cutout ^775- and ^850-band images of candidates. The image sizes are 2 arcseconds. 
The objects are named alphabetically (see Table 1) as marked at the upper-left corner of 
each panel of paired images. The images of the quasar, at the lower-right panel, are marked 
as "Q" . The band name is marked at the lower-left corner of each image. The contour at the 
center of each image marks the area where flux is collected for MAG_AUTO. The source color 
was calculated from MAG_ISO, which were derived from a compact core within the contour. 

Fig. 3. — Distribution of Z775 — ^850 color for simulated galaxy spectra at redshift z — 5.8. 
The open region represents LAEs, and the shaded region LBGs. Our selection of 1.3 < 
^775 — ^850 < 2.0 is expected to include most LBGs, but only about 1/3 of the LAEs. 

Fig. 4. — Field around object C in the ^850 band. Two additional red objects are marked. 
The angular scale of one arcsecond (~ 5.7 kpc of proper distance) is marked. 
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